Internal magnetic fields in thin ZnSe epilayers 
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Abstract 

Strain induced spin-splitting is observed and characterized using pump-probe Kerr rotation spec- 
troscopy in n-ZnSe epilayers grown on GaAs substrates. The spin-splitting energies are mapped 
out as a function of pump-probe separation, applied voltage, and temperature in a series of sam- 
ples of varying epilayer thicknesses and compressive strain arising from epilayer-substrate lattice 
mismatch. The strain is independently quantified using photo-luminescence and x-ray diffraction 
measurements. We observe that the magnitude of the spin-splitting increases with applied voltage 
and temperature, and is highly crystal direction dependent, vanishing along [110]. 
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ZnSe has been extensively investigated in the last decade for applications in optoelectronic 
devices operational in the blue spectral regiorT. More recently, the potential for use of ZnSe 
in spintronics devices has been evidenced by the long spin coherence times (~100 ns) in 
n-ZnSe and the ability to coherently transfer spins across a heterointerface between n-ZnSe 
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and GaAs."*" Experiments have further demonstrated electrical generation of spins in n- 
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ZnSe through both current-induced spin polarization and the spin Hall Effect at room 
temperature. 

Earlier studies in strained GaAs and InGaAs^ have revealed a A;-linear spin-splitting 
which acts on spins as an internally generated magnetic field . Despite the observation of 
electrically generated spins, no internal magnetic fields have been measured in ZnSe epilayers 
of micron-scale thickness. In this letter, we investigate the presence of electrically-induced 
internal magnetic field, Bi nt , in strained n-type ZnSe epilayers. We observe an effective Bi nt 
proportional to the drift velocity of electrons and note that, although the bulk inversion 
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asymmetry spin-orbit parameter in ZnSe is almost half that in GaAs, the spin-splitting 
energy scale is comparable. Coupled with the previous studies, this establishes the potential 
for developing all-electrical protocols for spin initialization and coherent manipulation in 
this wide band gap material. 

We measure three Cl-doped n-type ZnSe epilayers grown by molecular beam epitaxy on 
semi-insulating (001) GaAs substrates. The epilayer thicknesses are 100 nm, 150 nm and 300 
nm with carrier densities in the range of 4-5 x 10 17 cm~ 3 . Measurements of similar samples 
with 1.5-//m thick ZnSe epilayers and similar carrier concentrations show no evidence of 
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an internal magnetic fielcT . Figure la is an optical image of a patterned sample, showing 
channels etched along [110] and [110] directions. The channels have width w = 100 /im, 
length / = 235 fim along [110] (along [110], I = 350 jum), and thickness h = 100, 150 and 
300 nm. Annealed indium contact pads at the ends of each arm allow us to independently 
apply an electric field E along both orthogonal directions. 

Growth of ZnSe on GaAs causes in-plane compressive strain in the ZnSe film due to 
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the small lattice mismatch between the two zinc-blende semiconductors (0.25% at 300 K)r. 
X-ray diffraction (XRD) measurements have shown that although biaxial strain persists in 
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epilayers as thick as 1.6 /mf", the critical thickness beyond which dislocations begin to form 
in the ZnSe/GaAs system is about 0.17 /mi 10 . We quantify the strain in our samples using 
both photoluminescence (PL) and XRD measurements. Fig. lb shows PL at T = 50 K for 
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the 300-nm and 100-nm samples. Comparing the direct band gap of a 2-/mi ZnSe epilayer 
with similar carrier concentration at 50 K~ (dashed line in the figure) with the observed 
PL, we notice that the PL peaks are increasingly red-shifted relative to the unstrained 
band gap with decreasing epilayer thickness. By assuming that the PL peak occurs at the 
direct bang gap energy, we can use the shift in the PL to estimate the in-plane strain in 
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the samples - . To first order, these calculations predict compressive strain in the range of 
lxl0~ 3 to 3xl0~ 3 , which are comparable to the values quoted in existing measurements of 
samples of similar thicknesses 9 . XRD measurements (fig. lc and Id) show further evidence 
of monotonically increasing strain with decreasing epilayer thickness, again in agreement 
with previous studies. We do not observe any anisotropy in the strain about [115] and [115] 
directions in any of the three samples, and the diffraction pattern about the (004) peaks is 
found to be invariant with rotation about the growth axis [001], implying no tilt between 
the substrates and the epilayers. 

The samples are mounted in the variable temperature insert of a magneto-optical cryo- 
stat. The electron spin dynamics are probed with time-resolved Kerr rotation (KR) 
spectroscopy^ 3 "'^ in the Voigt geometry using pulses from a frequency-doubled mode-locked 
Ti:sapphire laser with a pulse duration ~150 fs and a repetition rate of 76 MHz. A circularly- 
polarized pump pulse normal to the sample surface injects electrons with spins polarized 
along the beam propagation direction z. After a time delay At, the Kerr rotation angle of 
a linearly polarized probe pulse measures the projection of the electron spin magnetization 
along z as it precesses about the transverse magnetic field B ext applied along y. The pump 
and probe energies are typically tuned between the range 2.70 - 2.80 eV, to maximize the 
KR effect, with powers of 1 mW and 100 /xW, respectively. The beams are focused to a 
spot size of ~20 /im, and the circular polarization of the pump beam is modulated with a 
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photoelastic modulator at 50 kHz for lock-in detection. Time-resolved measurement 3" have 
shown that the electron g-factor for ZnSe in these samples is 1.1. 

In order to characterize the internal magnetic field Bi nt , we measure KR as a function 
of the applied magnetic field B ext at a fixed pump-probe delay time of At = 13 ns while 
varying the applied electric field E and pump-probe separation d. Fig. 2a shows the data 
from the 100 nm sample (with the largest strain) at T = 125 K when the pump and probe 
are spatially overlapped (d = 0). The signal is an oscillatory function of the total magnetic 
field B = B ext + Bi n t and is given by Acos(gfiBBAt/hp^ where A is the amplitude, g the 
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effective g-factor of the sample, [lb the Bohr magneton and B the magnitude of the total 
magnetic field experienced by the electrons. For KR along [110] in the absence of an applied 
electric field (top trace), the average k=0 and therefore B int = and the oscillations are 
centered at B = B ext = 0. For E parallel to B ext , the center peak is suppressed (center 
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trace). In this geometry, B int is perpendicular to B ext ~* resulting in a total field magnitude 
of B = \J B\ xt + Bf nt which is always greater than zero for a non-zero Bi nt . For larger 
applied voltage, B int increases, resulting in greater suppression of the central peak. When 
the sample is rotated so that E is perpendicular to B ext (bottom trace), both the external 
and internal magnetic fields are along y and add directly so that B = B ext + Bi nt and the 
oscillatory signal is centered at -B int . It is noted here that all of these effects are observed in 
the channel along [110] direction only. For the channel along [110], surprisingly no internal 
fields in any geometry are observed in any of the samples. 

In both cases above, with non-zero electric field the amplitude of the signal decreases 
with increasing voltage, which is further investigated by spatially separating the pump and 
the probe by a distance d along the direction of E (Fig. 2b). Due to the laser profile of the 
pump beam, the optically injected spins have a Gaussian spatial profile which is centered at 
d = fj>m when E — 0. An applied voltage (E ^ 0) imparts a non-zero average momentum 
k to the injected spin packet, causing it to drift with an average velocity v s . We fit the 
amplitude A as a function of d to determine the center position d c = 6 /mi of the spin 
packet with applied voltage of 20 V after 13 ns, giving v s = 0.46 /xm/ns. Fig 2c shows the 
spatial variation of B int throughout the spin packet for E parallel to B ext . Spins at the 
leading edge of the packet experience a larger Bi nt than the trailing edge. This variation is 
due to the spread in the drift velocities of the spin packet arising from spin diffusion. The 
reported value of B int for each voltage are obtained from a linear fit at d — d c . 

Fig. 3a shows both data and fits for A (top) and B int (bottom) of the spin packet as 
measured by KR at different applied voltages in the 100 nm sample at T = 50 K. Bi nt 
increases linearly with applied voltage (inset, top) confirming the ^-linear dependence. It 
also increases with T, almost doubling as T changes from 25 to 150 K (inset, bottom). It has 
been observed^ that the in-plane biaxial strain in a 100 nm ZnSe/GaAs epilayer increases in 
magnitude with increasing T due to the thermal changes in the lattice and elastic constants 
of both epilayer and substrate, which may explain the changes observed here. However, 
without more systematic data on the alteration of Bi nt with strain, this remains speculative. 
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While spin precession in these samples persist to higher temperatures, at T > 150 K the 
spin lifetime is not long enough to quantify the internal field. 

Fig. 3b shows the variation of the spin-splitting energy A = g^iB^mt with the spin drift 
velocity v s for all samples at T = 50 K. From linear fits to the data, we use the slope j3 = 
A/v s to characterize the samples and find (3 = 101, 86 and 103 neV ns//im for the 300, 150 
and 100 nm samples, respectively. Previous measurements^ in GaAs samples with strain in 
a similar region of 2-3xl0 -3 have shown (3 values ranging between 75-125 neV ns/ /jhi. 

The contact resistances in the samples are non-negligible and the effective electric field 
across each channel cannot be simply calculated from the applied voltage. Resistivity ob- 
tained from transport measurements is used to calculate the electric field in the two thicker 
samples, and in Fig. 3c we compare B int between the 300 and 150 nm samples as a function 
of E. Surprisingly, we observe the thicker sample, with less strain, having a larger internal 
magnetic field for the same electric field. Unfortunately, for the 100 nm sample the contact 
resistance is too large to allow any reliable estimate of the electric field 

In conclusion, we have quantitatively measured the ^-linear internal magnetic field as a 
function of voltage in three samples. At this point we do not have a complete understanding 
of the correlation between the internal magnetic field and strain in ZnSe since comparison 
between the 300 and 150 nm samples show an increased spin splitting with decreasing strain. 
This effect is not systematic, as we do not observe any internal fields in the largely strain 
relaxed 1.5 /xm thick sample. We also observe Bi nt only along the [110] direction in the 
samples while the strain is isotropic, which, though not well understood, may be exploited 
to design spin manipulation using topological effects 1 ^ . 
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FIG. 1: (a) Optical image of processed sample showing the crystal orientation of the epilayer 
(/ = 235 /im, w = 100 /im). The coordinate axes define the measurement geometry, (b) PL 
for 300 and 100 nm samples showing the strain-induced shift. Arrows indicate peak emission 
energies. Dashed line indicates the bulk, unstrained energy gap. XRD spectra for the (c) (004) 
and (d) (115) reflections. The appearance of interference fringes in the (004) reflection in the 100 
nm sample indicate phase coherent growth and may be used to further verify the epilayer thickness. 
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FIG. 2: Data from 100 nm sample at T = 125 K. (a) Kerr rotation with pump-probe overlapped 
(d=0) at At = 13 ns with (top) E = 0, (center) E \\ B ex t and (bottom) E _!_ B ex t- (b) Schematic 
of measurement geometry for spatially resolved scans, (c) A and Bi n t, obtained from fits described 
in the text, as a function of pump-probe separation d with E\\ B ex t- 
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FIG. 3: (a) A and Bi n t (top and bottom) as a function of pump-probe separation d with E\\B ex t, for 
a series of applied voltages V at T = 50 K in 100 nm sample. (Inset, top) Bi n t varies linearly with 
applied voltage. (Inset, bottom) Bi n t as a function of temperature T for V = 20 V. Comparison 
of (b) electric field induced spin-splitting A at 50 K as a function of the spin drift velocity v s for 
all samples and (c) Bi n t for 300 and 150 nm samples varying with E at 50 K. 
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